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Abstract 

This paper develops approximation methods for path-dependent functionals, which have been 
used in many applications involving path-dependent objective functions. By a first glance, the 
problem may appear as a standard approximation of a stopping time problem using traditional 
techniques. Nevertheless, a closer scrutiny would reveal that the problem under consideration 
is far more challenging and difficult. The difficulty is mainly because the traditional dynamic 
programming approach falls apart, not to mention any hope for a closed- form solution or any 
viable numerical PDE methods. This paper focuses on approximation methods in Monte Carlo 
simulations. To the best of our knowledge, convergence analysis is available only in few special 
cases due to the complexity of the nature of path dependence. In contrast to the traditional 
approach, this work provides a non-traditional convergence method in Monte Carlo analysis, 
and establishes a general framework for the convergence analysis using the Markov chain ap- 
proximations. In particular, the approach is based on actual computations under the Skorohod 
topology. Some examples such as the approximation of discretely monitoring barrier option are 
considered. 
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1 Introduction and Examples 



In many applications, one needs to evaluate path-dependent objective functions. Many of such cal- 
culations require the handling of stopping times in which the stopping boundaries are non-constant 
but are path dependent. The applications stem from evaluations of functionals in networked systems 
as well as in finance. In addition, they also arise in Euler's approximation to solution of stochastic 
differential equations. This paper is concerned with approximation methods for computing such 
objective functions. By a first glance, the problem may appear as a standard approximation of 
a stopping time problem involving traditional techniques. Nevertheless, a closer scrutiny would 
reveal that the problem under consideration is far more challenging and difficult. Even though the 
problem may still be viewed as stopping time problems, the targets for the stopping times are not 
fixed constants as most of those treated in the literature, but rather a path-dependent criteria. 
This is where the main difficult arise because the traditional dynamic programming and Markovian 
properties fall apart. 

To begin, there are virtually no closed-form solutions for the path-dependent objective functionals. 
One has to look for alternatives. The next possibility is numerical approximation using partial 
differential equation (PDE) based techniques. Nevertheless, such possibility is ruled out due to 
the lack of Markovian properties. The only often used technique left is the Monte Carlo method. 
While most of the existing methods for treating Monte Carlo procedures are somewhat ad hoc, this 
work develops a systematic alternative method for analyzing the convergence of the approximation 
algorithm. In this paper, we use a weak convergence approach and carry out the convergence 
analysis under the Skorohod topology. One of the main ingredients is the use of a generalized 
projection operator. With the use of such projections, we proceed to by exploring the intrinsic 
properties of the Skorohod space. The convergence analysis developed in this paper provides a 
thorough understanding of the nature of path dependence and a general framework for handling 
many path-dependent problems. 

1.1 Path-dependent Objective Functions 

We are interested in approximating path-dependent functions in a finite-time horizon. For sim- 
plicity, we focus our attention on the time interval [0, 1] to avoid using more complex notation. 
Let C[0, 1] be the collection of continuous real-valued functions defined on [0,1], and D[0, 1] be 
the collection of all RCLL (right continuous with left-hand limits) functions on [0, 1] . With 
F = {Tt : t G [0,1]} and J- = J-i, let F, F) be a filtered probability space satisfying the 

usual conditions; see [10]. As a convention, we use capital letters to denote random elements and 
use lower case letters for deterministic (or non-random) elements. For instance, x G -^[0, 1] is a 
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RCLL path defined on [0, 1], while X : i— )■ D[0, 1] is a random RCLL process with sample paths 
in D[0, 1]. As usual, X can be regarded as a function of two variables time t and sample point 
uj. That is, for each fixed u £ X{-,uj) denotes a sample path and for each fixed t, X{t, •) is a 
random element in R. 

Throughout the paper, a,/3 G C[0, 1] are given satisfying a{t) < I3{t) for all t G [0, 1]. The domain 
of our interest is defined by all the points bounded by a and /3, that is, 

Q = {{a,t): a{t) <a< f3{t),t G [0, 1]}. 

The cross-section Q{t) = {a{t), P{t)) C M may be time-dependent. We also denote its boundary by 

dQ = {{a{t),t) : t G [0, 1]} U {(/3(t), t) : t G [0, 1]}. 

We study an F-adapted continuous process X : [0, 1] x 17 i— t- C[0, 1]. Let r be the first hitting time 
to the boundary of the time-dependence domain 

T = inf{t > : X{t) ^ Q{t)} A 1. (1) 

Given a positive integer m, let (ci, C2, . . . , c^)' G [0, 1]*" be a column vector, representing a vector 
of scaling factors. Define a set of random times 

T = {Tj = CjT : 1 < J < m}. 

Note that random times Tj are measurable with respect to but may not be a stopping time 
whenever < Cj < 1. Now we introduce the notion of generalized projection operator by defining 
n as 

UriX) = {X{Tj) : 1 < i < m} G M"". 

The n is a generalized projection in the sense that it distinguishes from the conventional projection 
of a path by projecting X to a collection of random times instead of fixed times; see [3]. 

For X G D[0, 1], we define a maximum process x* by x*{t) = supo<<(<t It is easy to check that 
X* G -D[0, 1]. For some measurable function g : ^ M, we are interested in the computation of 

the objective functional V : C defined by 

V = E[g{Ur{X),X*{T),T)]. (2) 

The value function may be a measurable function of the initial states V{x, t) if the functional V is 
given as a non-path-dependent form of supE[/(X(T)] for some function / : R i— )• M. However, we 
emphasize that our main interest is the computation when the functional is path-dependent in the 
form of ([2]), which has a wide range of applications in mathematical finance. 
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Remark 1 Related literature in connection with a hitting time under a Markovian framework 
can be found in Dufour and Piunovskiy [6], de Saporta, et al. [5], and Szpruch and Higham 
|18j among others. In particular, in [6], the existence of optimal stopping with constraints is 
established using a convex analytic approach. A numerical method for optimal stopping of piecewise 
deterministic Markov processes is considered in [5]. Bounds for the convergence rate of their 
algorithm are obtained by introducing quantization of the post jump location and path-adapted 
time discretization grids. In [18], an application in physics (thermodynamic limit) involving mean 
hitting time behavior is considered. 

1.2 Examples 

The evaluation of path-dependent objective functions is of great interest in many networked systems 
as well as in financial applications such as the option pricing. To illustrate, consider the underlying 
stock price X that is a non-negative martingale process under P that is the corresponding equivalent 
local martingale measure. Assuming zero interest rate, then the value y of ([2]) can be considered as 
a general form of a class of option prices including look-back option, rebate/barrier option, Asian 
option, Bermuda option etc. For illustration purpose, let us fix 

T = I — T : 1 < i < ml. 
I m J 

1. (Barrier option) Let —a{t) = oo and /3(i) = 1, and thus r = inf{s > t : X{s) ^ (— oo, 1)} A 1. 
The price formula of the up-and-in barrier call with strike 1/2 at maturity T = 1 is given by 

v = n{x{i)-\)^hr<i}\- 

In this case, the payoff function is 
which is discontinuous. 

2. (Discretely monitoring barrier option) As mentioned in [U and [12], treating a class of so- 
called monitoring barrier options, most models assume continuous monitoring of the barrier 
like the above. However in practice most, if not all, barrier options traded in markets are 
discretely monitored. Unlike their continuous counterparts, there is essentially no closed form 
solution available, and even numerical pricing is difficult. For instance, again we set the same 
barriers —a{t) = oo and /3(t) = 1. However, we will replace r in the above continuous time 
barrier option by 

f = inf{s = i/m : X{s) ^ (-oo, 1), i = 1, 2, . . . m} A 1, 
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The price formula of up-and-in barrier call with strike i is then 



V = E 



X(l)-l] I. 



2J -'{^<i}J- 

Observe that V can be rewritten as V = E[(X(1) — ^)'''/{max(n7-(X))>i}]) and the payoff 
function g is accordingly 

/ 1\ + 

g{xi, X2, ■ ■ ■ , XjYi, X,iji^l,t) = \^^~2J ^{max{a;i,...,a;m}>l} • 

Note that g has linear growth, is unbounded, and is discontinuous. 
1.3 Computational Methods 

To evaluate the path-dependent objective functions, a major difficulty arises from the lack of 
Markovian conditions. This rules out the possibility of any PDE based methods including the 
finite difference or finite element methods. In addition, the time-dependent barriers a and /3 also 
create additional layer of difficulty. In this paper, we consider Monte Carlo methods to obtain 
a feasible estimation of the path-dependent problems. Monte Carlo method is a class of compu- 
tational algorithms that relies on repeated random sampling to evaluate its deterministic value 
using its probabilistic fact. This includes Euler-Maruyama approximation [11] and Markov chain 
approximation [8], [13], and [TJ] among others. 

The general idea of the Monte Carlo method in this vein is the following. For each sample point 
w G il, use X'^{-,uj) G D[0, 1] to denote a simulated path for the underlying process X{-,u) G C[0, 1] 
by a certain Monte Carlo method with a small parameter h. Define the approximated stopping 
time and projection operator by 

= inf{t > : X^{t) i (a(t),/3(t))} A 1 and = {t^ = cjt^ : 1 < i < m}. 

One can approximate V in this way by computing 

= E[g{Ur4X%X^^*{T%T% (4) 

where 

Urn{X^) = {X\T!;):l<j<m}. 

The goal is to design an appropriate Monte Carlo method so that the desired convergence takes 
place eventually, i.e., 

Hm V'' = V. 

As the objective value is given as an expectation form of ([2]), V is invariant under the same 
distribution, and the usual requirements for the Monte Carlo method are intuitively given as follows: 
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(HI) X converges in distribution to X as /i — )• 0, denoted by X =^> X. 
(H2) g is continuous. 

A couple of natural questions are as follows. 

• Are (H1)-(H2) sufficient to guarantee the desired convergence lim/j_).o V^'* = V? 

• Can (H2) be possibly weakened to some discontinuous function g so that the barrier option 
pricing (see Section II. 2p can be included? 

1.4 Counter Examples 

Interestingly, there exist circumstances that lead to counter examples in connection with the desired 
convergence under (H1)-(H2). 



Figure 1: Demonstration of tangency problem: The sample paths of the two processes (solid line 
and dotted line) are fairly close to each other, but their first exit times are far apart indicating a 
lost of continuity. 

A noticeable counter example is given in Figure [H which is motivated by the so-called tangency 
problem. Consider two underlying processes X^ (solid line) and (dotted line) in Figured) No 
matter how close X^ and X'^ are at the initial states, the difference between their first exit time ri 
and T2 could be very far away; see [2]. The above idea is formulated as the following example. 

Example 1 Let {X{s) : < s < 1} be a deterministic process given by X{s) = 1 - (s - i)2. Let 
a{s) = — oo and /3s = 1 for all s > 0. Then the exit time of X is 



o 




T = inf{s > : X{s) ^ {a{s), f3{s))} A 1 = ^. 
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Define a family of the processes parameterized by h with X (s) = X{s) — h. Although X converge 
to X in L°° as /i — )■ 0^ , we note that 

t'' = inf{s > : X'^(s) ^ (a(s), /3(s))} A 1 = 1 

not converging to r = 1/2. 

The next example is different from the first one in the sense that there is no obstacles and the 
underlying process is taken from a Bessel process. This system is often used in finance to model the 
dynamics of asset prices, of the spot rate and of the stochastic volatility, or as a computational tool. 
In particular, computations for the celebrated Cox-IngersoU-Ross (CIR) and Constant Elasticity 
Variance (CEV) models can be carried out using Bessel processes. 

Example 2 Let X{t) be the Bessel process of dimension 3 with initial X{0) = 1, i.e., 

X«) = l + «'«)+/'^, (5) 

which is a well known strict local martingale with E[X(1)] < 1; see Page 336 of [9]. Define 
X^{t) = X{t) A Fix an arbitrary partition Tn = {ti < ■ ■ ■ < tn} C [0,1], and an open set 
O C [0, iV]". Note that, for any h £ {0,1/N), we have 

p{nr„ {x^) G 0} = P{nr„ (x) GO} - p{nr„ (x) G o, nr„ ixn>l} 

and 

<P{nr„(X)GO,nrjX*)>i} 
<P{nr„(X*)> i}^0 as/i^O. 

This implies limh^oF{Ur^{X'') G O} = F{Ur„{X) G O}, and thus Ur„{X'') Ur„{X) due 
to the arbitrariness of the set O. Moreover, Kolmogorov consistency theorem together with the 
uniqueness of the weak solution of SDE ([5]) implies that X^ =^ X due to the arbitrariness of n and 
Tn,. However, since X^ is a bounded local martingale, hence a martingale, we have the following 
undesired convergence: lim/j^o IE[X'*(1)] = 1 > E[X(1)]. 

1.5 Goal and Outline of the Paper 

The above examples show that (H1)-(H2) may not be sufficient for Monte Carlo simulation to be 
convergent to the right value of ([2]). This propose the following question: 

(Ql) Given X^ =^ X, what is a sufficient condition for the convergence lim.h-).oV^ = VI Is it 
possible to weaken the continuity of g to cover the barrier option? 
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In this paper, we first explore a rigorous proof of the convergence in Section [2j The results in 
Theorem [6] provide a set of sufficient conditions for the convergence. It can also explain why 
Examples [J and [2] do not work. 

Our approach is based on actual computations using the Skorohod metric in the D[0, 1] space. This 
is the first attempt in this context to the best of our knowledge. Such an approach is advantageous. 
For example, different types of options mentioned the convergence of Monte Carlo methods have 
been studied separately; see [5] and the references therein. However, Theorem [H] provides a unified 
theoretical basis for the convergence in a general path-dependent form with possibly discontinuous 
payoff function g and finite time-dependent barrier This covers a wider range of 

applications of Theorem [H] to various types of options, for instance, barrier options, Asian options, 
look back options, etc. It is also notable that the assumptions on the barrier in the path-dependent 
problem has some similarities to the assumptions imposed for stochastic control exit problems; see 
[2] . To complete the approximation of the value F of ([2]) , one should consider the following question 
in addition to (Ql), 

(Q2) Given X, how does one construct approximating processes X'^ such that X'^ =^ Xl 

Among many of the possible answers to (Q2), we will mainly focus on the construction of X^ in 
the framework of Markov chain approximation in Section [3j Recall that from the book |14j , a 
family of continuous approximating processes X^ constructed using Markov chain is convergent to 
X in distribution if the Markov chain satisfies the local consistency [141 equation (9.4.2)]; see also 
\14:\ Theorem 10.4.1]. Compared to [13], we investigate the weak convergence for more generalized 
local consistency condition. This enables us to cover various types of Monte Carlo simulations in 
the framework of Markov chain approximation to verify its weak convergence using the generalized 
local consistency. It essentially extends the use of Markov chain approximation. 

To proceed, the rest of the paper is arranged as follows. Section [2] provides sufficient conditions 
for convergence of our procedure. Then the convergence of the algorithms is established. Section 
E] presents the weak convergence of Markov chain approximation under the generalized local con- 
sistency. Also provided here is the answer: Can we get strong convergence? Finally, the paper is 
concluded with some further remarks in Section HI 
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2 Sufficient Conditions for Convergence 



2.1 Preliminaries 

We use the notation given in [3]. Define a metric on D[0, 1] by 

\\x-y\\= sup \x{t) -y{t)\, \/x,y e D[0,l\. 

te[o,i] 

Then, the continuous function space on [0, 1], denoted by C[0, 1], is complete with respect to the 
uniform topology with the above metric || • ||. On the other hand, the RCLL function space defined 
on [0, 1], denoted by -D[0, 1], is equipped with the Skorohod topology with the metric 

\\x - y\\s = inf {||A - /||, ||xo A - y\\}, \/x,y G D[0, 1], 
aga 

where x o y denotes the composite function of x and y, and A is the collection of all continuous 
increasing functions A on [0, 1] with A(0) = and A(l) = 1. Note that D[0, 1] is not complete under 
the metric || • but it is complete under an equivalent metric || • ||° defined by 

m - m 



inf (sup log^^^ ^ yx,yeD[0,l]. 

XeA L s<t t — s ) 



\\x-y\\, ^ ^ . . 

AeA L s<t 

For notational convenience, we will adopt || • \\s in the rest of the paper. Recall that, for x G -D[0, 1], 
we define x* G D[0, 1] by x*{t) = supo<s<t 

Lemma 2 Let vr* : I? [0,1] i— )■ I? [0,1] hy 7r*(x) = x* . Then, n* is continuous with respect to the 
Skorohod topology. 

Proof: Since A is strictly monotone, x* o A = (x o A)*. Therefore, 
\\x*-y*\\s = inf {||A-/||, sup |x*(A(t)) - 

tG[0,l] 

= inf{||A-/||, sup \{xoX)*{t)-y*{t)\} 
te[o,i] 

= inf {||A - /||, sup ||xo A-y||[o,t]} 
te[o,i] 

= ~ -^11' Ik o ^ - 2/11} = Ik - y\\s- 

aga 

So, vr* is continuous in D[0, 1] with respect to the Skorohod topology. □ 

2.2 Main Convergence Results 

We make the following assumptions. 

(Al) If /? < oo, then F{iX - /3)*{1) > 0} > 0. If a > -oo, then F{{-X + a)*(l) > 0} > 0. 
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Note that Condition (Al) means that the process hits both barriers a and /3 with positive prob- 
abihty, whenever — a or /3 are bounded above. This is not a restriction, since one can simply set 
a{t) = —CO (resp., /3(t) = oo) for all t G [0, 1], if X never hits a (resp., /?) almost surely in P. If 
—a{t) = f5{t) = oo, then r = 1. 

(A2) X satisfies inf{t > r : X{t) ^ [a(t), = r almost surely P. 

Condition (A2) requires that the boundary dQ is regular with respect to the process X. Note that, 
for any small e > 0, X exits from Q in the interval (r, r + e) under (A2). Loosely speaking, the 
condition (A2) means that the process X{t) exits Q immediately after it hits the boundary at r. 
More discussions are referred to Remark [71 
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Figure 2: Illustration for the partition of C[0, 1] 

To proceed, let's partition the space C[0, 1] as follows: Define 

Ci = {x E C[0, 1] : T^[x) < 1, x{ir{x)) = ^iirix)), M{t > ir{x) : x{t) > p{t)} = tt{x)}, 
C2 = {x e C[0, 1] : 7r(x) < 1, x{tt{x)) = a{iT{x)), inf{t > tt{x) : x{t) < a{t)} = vr(x)}, 

and 

C3 = {x G C[0, 1] : 7t{x) = 1}, and C4 = C[0, 1] \ (uf^iQ). 

Accordingly, we can write C[0, 1] = U^^^^Cj and Cj H Cj = for i / j. As for the illustration, one 
can see that the four curves depicted in Figure [2] belong to four different subsets separately, that 
is, li G Ci for i = 1,2,3,4. Note that (A2) implies that 

F{X e ^f=iCi} = 1 or equivalents F{X ^ C4} = 0. 
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Lemma 3 Suppose (Al) and (A2) hold, and X. Then, {X^,t^) (X,r) as h^O. 

Proof: Let vr : -D[0, 1] ^ [0,1] be the mapping 7r(x) = infji : x{t) ^ /3(t))}, and vr : 

D[0, 1] —7- D[0, 1] X [0, 1] defined by 7r(x) = {x,Tr{x)). Note that we want to show that, under (Al) 
and (A2) 

X''^ X hnphes that fr{X'') Tr{X). 

Since (A2) imphes that F{X G uf^j^Cj} = 1, by the continuous mapping theorem [3l Theorem 2.7], 
the conclusion holds provided that the mapping vr is almost surely continuous under the probability 
induced by FX~^, or equivalently vr is continuous at arbitrary x G uf^j^Cj. Hence, it's enough to 
prove that vr is continuous at arbitrary given x G uf^^Ci with respect to the Skorohod topology in 
the rest of the proof. 

We assume —a{t) + f3{t) < oo for each t without loss of generality. Fix x G uf^j^Cj, and be 
an arbitrary sequence in D[0, 1] satisfying — x\\s — >■ as n — )• oo. By the definition of || • \\s, it 
implies that there exists A„ G A such that 

||A„ — /|| — )■ 0, \\xn o A„ — x\\ — )• 0, as n — )• oo. (6) 

• (Case 1). Suppose a; G Ci, and define, for any e > 

5e := I sup {x{t) - /3(t)). 

^ 7r(a;)<t<7r(a;)+e 

By the definition of Ci, we have 5^ > for all e > 0. Moreover, by continuity of x, there 
exists a time G (vr(x), Tr{x) + e) such that x(t£) > /3(te) + 6s- Next, we observe that 

\xniKite)) - x{ts)\ < \\xn o A„ - x|| -^0, as n OO. 
This implies that for some large Ns^i 

Xn{Xn{te)) > x{ts) - y > /3(t,) + ^, Vn > iV,,i. (7) 

On the other hand, there exists > due to the continuity of /3 such that 

sup |/3(t)-/3(t,)|<^. 

Moreover, since ||A„ — /|| — s- 0, we have for some large enough Ns^2, 

\\\n - I\\ < re, Vn > Ns^2- 
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Therefore, it leads to 

|/3(A„(t,))-/3(t,)| < sup |/3(t)-/3(t,)| Vn>iV,,2. (8) 

([7]) together with ^ imphes 

Xn{K{te)) > PiKite)) + ^, Vn > iV, := max{A^,,i , iV.^a}, 

or equivalently, 

vr(x„) < Xn{te) < A„(7r(rE) + e), Vn > N^. 

Finally, taking the lim sup„_j.o(3 on each side of the above inequality, and using \\Xn — /|| — > 
of (HI), we have 

lim sup 7r{xn) < lim sup A„(7r(2;) + e) = 7r(x) +e. 

n—^oo n—^oo 

So we conclude that by arbitrariness of e 

lim sup vr(x„) < tt{x). 

n— >-oo 

• (Case 2). Again suppose x £ Ci, and we prove the reverse inequality. By the definition of 
the first hitting time, for any e > 0, there exists 6^ > such that 

- x)*(7r(x) - e) > 6e, (x - a)*(7r(x) - e) > 5,. (9) 

Also, dll) implies that, \\a o A„ — a\\ — )• 0, o A„ — /3|| — )■ and o A„ — x\\ — )• 0. Using the 
triangle inequality, we have 

\\{xn - a) o A„ - (x - a)\\ ^ 0, ||(x„ - /?) o A„ - (x - /3)|| ^ 0. 

So, there exists a large such that 

||(x„ - q) o A„ - (x - a)\\ < |, ||(xn - /3) o A„ - (x - /3)|| < Vn > iV,. 
In particular, this leads to 

S 6 

{a - Xn) O Xnit) < (a - X)(t) + ^, {Xn " /3) o Xn{t) < (x - /3)(t) + y , G [0, 1], 71 > N,. 

By taking supp j^^^^.).^] on both sides of the above inequalities and using ([9]), we have 

S S 

{xn - /3)*A„(7r(x) -e)<{x- /3)*(7r(x) - e) + ^ < Vn > iV„ 

and 

(q - Xn)*A„(7r(x) - e) < (q - x)*(7r(x) - e) + ^ < -y, Vn > TV^. 
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Thus, we have 

7r(x„,) > A„(7r(2;) -e), Vn > N^. 
Taking hminf„ and using ([6|) and the fact of the arbitrariness of e, we have 

hminf7r(x„) > 7r(x). 

n 

Summarizing the above, we have proved the continuity of vr at x G Ci with respect to the Skorohod 
topology, i.e., hm„_s.oo vr(x„) = 7r{x) whenever x„ — )• x S Ci with respect to the Skorohod topology. 
Similar arguments as in Cases 1 and 2 yield that vr is also continuous at x G C2. Moreover, vr is 
continuous at x G C3 thanks to Lemma [2l This completes the proof. □ 

Lemma 4 Suppose (Al) and (A2) hold. If{X^,T^) {X, t) as h ^ 0, then 

1. {X^,X^'*,T^) {X, X*,t) ash^ 0, 

2. iU^hiX''),X'''*{T''),T^) ^ {UriX),X*iT),T) ash^O. 

Proof: By Lemma [3 the mapping ttI : D[0,1] x [0,1] ^ -D^[0, 1] x M by 7rl{x,t) = {x,x*,t) 
is continuous in the product space in the Skorohod topology, where £'^[0,1] = £'[0,1] x D[0,l]. 
Therefore, (X'^,r'^) ^ (X,r) implies 7r|(X^,T^) ttI{X,t) by mapping theorem, which in turn 
concludes the first result {X'' , X'''* ,t'') =^ {X,X\t). 

Next, for simplicity, we only consider the special case with (01,02) = (1/2,1), so that the desired 
result is reduced to show that 

{X^{t''),X'^'*{t''),X\t^/2),t^) {X{t),X*{t),X{t/2),t), ash^O. 

Let 7f : D^[0,1] x [0,1] ^ by 

7t{x,y,t) = {x{t),y{t),x{t/2)),y{x,y,t) G D^[0,1] x [0,1]. 

Fix {x,y,t) such that (x, y) is continuous at t and t/2, respectively. Suppose {xn,yn,tn) {x,y,t) 
in their product space. Then, there exists A.„ G A such that 

||A„ - /|| 0, ||xn - X o A„|| + ||?/„ - yA„|| + |t„ - t| 0, as n 00. (10) 

Therefore, 

\Xn{tn) - x{t)\ < |x„(t„) - x(A„(tn))| + |x(A„(t„)) - x{t)\ 

< \\Xn - X o A„)|| + |x(A„(t„)) - x{t)\ 
Note that the first term ||xn — x o An|| — t- as n — )• 00 by (jlOp . On the other hand, since tn ^ t 
and ||A„ — /|| — )• as n — )• 00, we have |A„(t.„) — t\ — )• as n — )■ 00. Hence, the second term 
\x{Xn{tn)) — x{t)\ — )• by continuity of x at t. So we have Xn{tn) — ^ x{t). Similarly we can show 
Vnitn) — ^ uit) and x„(t„/2) — )• x{t/2). Thus, we conclude 
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vr is continuous at {x,y,t) whenever {x,y) is continuous at t and t/2. 



Note that, {X,X*) is continuous on [0,1] almost surely in P. In particular, X{-,u}) and X*{-,uj) 
are continuous at t{u) and t{uj)/2 almost surely in uj with respect to P. Therefore, if we set 

A := {(jj ^ Q. : {X{-^uj),X* {-^oj)) is continuous at t{uj) and r(a;)/2}, 

then it's equivalent to P(j4) = 1. Therefore, with random element Z : 17 i— )■ (X, X*,r), we have 
PZ~^(_D^[0, 1] X [0, 1]) = ¥{A) = 1, and vf is continuous almost surely in PZ~^. This also implies 
the mapping tt : L»2[o, 1] X [0,1] ^ by 

7r(x,y,t) := x(t/2), t) = (7r(x, y, t), t), V(x, y, t) e Z)2[0, 1] x [0,1] 

is continuous almost surely in ¥Z~^. Applying mapping theorem [3l Theorem 2.7], together with 
the previous result of {X^ , X^'* ,t^) =^ {X,X*,t), we conclude the second result. □ 

Next, we make an assumption on the payoff function g. 

(A3) g : M'""''^ — )■ M is an almost surely continuous function with respect to PZ~-^, where Z := 
(H-j-iX), X* {t),t). Moreover, g satisfies either one of the following growth conditions: 

1. g is a bounded function; 

2. g is a function with linear growth and {X^(t) : h > 0,t £ [0, 1]} is uniformly integrable. 
The following result is a consequence of Lemma HI 

Proposition 5 Assume (Al)-(A3). Let X be a Tt-adapted continuous process with initial X[{)) = 
X, and X^ he a sequence of RCLL processes satisfying X^ ^ X as h ^ 0. Then, lim/i_>.o = V. 

Proof: Since (Al) and (A2) hold. Lemma H] implies that 

iUr4X^),X^'*{r^),r^) iUriX),X*iT),r) 

as /i 0. Together with (A3), it results in liuih-^oV^ = V; see O P. 25, 31]. □ 
Now, we are ready to answer question (Ql). 

Theorem 6 Suppose (A1)-(A3) hold. If X^ ^ X ash-^0, then lim/,^o V''' = V. 
Proof: This is a consequence of Lemmas O HI and Proposition [H □ 

Theorem [6] holds under assumptions of (A1)-(A3). Recall that (Al) is not a restriction. We will 
discuss (A2), (A3), and the uniform integrability in what follows. 
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Remark 7 (Discussions on (A2) and Example [T]) In fact, (A2) is a requirement on the reg- 
ularity of the boundary dQ with respect to the process X, and it is referred to as r'-regularity for 
simplicity; see [T7j. Note that since X in Example [1] violates r'-regularity (A2), by observing 

inf{t > T : X{t) i [a{t),P{t)]} A 1 = 1 > 1/2 = r, 

it yields the convergence to the wrong value. In other words, (A2) is crucial for the investigation 
of the convergence. 

Remark 8 (Discussions on (A3) and Example [2]) Assumption (A3) is the requirement on 
the function g. First of all, it allows discontinuity of but it can not be too much discontin- 
uous in the sense that it is at least required to be almost surely continuous. However, it is already 
enough to include option pricing for the discontinuous payoff, such as barrier option. 

In (A3), another issue yet mentioned is the growth condition of g. In particular, if g is of linearly 
growth function of the underlying price like in the call type option, then one shall verify uniform 
integrability. For instance. Example [2] converges to a wrong value since {X{t) : t G (0,1)} is not 
uniformly integrable while the payoff function g[x) = x is linear growth. Recall the definition of 
uniform integrability. A set of random variables {Ky : 7 E F} is said uniformly integrable, if for 
any e > 0, there exists a compact set K such that 

supE[y^/{y^^^}] < e. 
7er 

Next, we present a proposition to ease the verification of the uniform integrability, which will be 
used later practically. The proof is referred to [7j. 

Proposition 9 A family of real-valued random variable : 7 G F} for some index set F is 
uniformly integrable z/ sup^ E[|Ky|P] < 00 for some p> 1. 

3 Weak Convergence of Markov Chain Approximation 
3.1 Markov Chain Approximation 

In this section, we establish the weak convergence of Markov chain approximation for multidimen- 
sional stochastic differential equations. Let Y = {Y{t) : t G [0,1]} be the unique weak solution 
of 

dY{t) = b{Y{t),t)dt + a{Y{t),t)dW{t)- y(0) = y, (11) 
where b : R'^+^ ^ R'^, and W is a standard M'^i Brownian motion, and a : R'^^^ 1— M'^^'^i. 



15 



Let tQ = < ti < ■ ■ ■ < t'^ = 1 he a sequence of increasing stopping times with respect to a discrete 
filtration {J"/^ : i = 0, 1, . . .}, and : i = 1,2, . . . N} he a sequence J-j^ measurable Markov chain 
in R'^ with transition probability 

F{Y.!l, G dylY/' = X, = t] = p^{t, x, y). 

We use = {Y^{t) : t G [0, 1]} to denote piecewise constant interpolation 

y\t) = Y.^l'ht-<t<t'i,,y (12) 

For notational simplicity, we set At^ = t^j^^ — and AY^ = Y^^^ — Y^. The interpolation of the 
Markov chain process Y^ is said to be locally consistent, if 

(LCI) E[sup„AtJj] = 0{h), 

(LC2) mY.!:\Y^^ = = i] = E[At::|y„'^ = yA = t\- {b{y,t) + 0{h)), 

(LC3) cov(Ay„'^|K„'^ = = t) = E[Atf5|y„^ = yA = t]- + o{h)). 

In local consistency (LC2) and (LC3), 0{h) is understood as a d-dimensional vector or d x d 
dimensional matrix with each element being 0(h), whenever it is appropriate. 

(A4) b and a are Lipschitz in y and Holder- 1/2 continuous in t, i.e., with <j) = b, a, 

\^{yi,ti) - 0(2/2,t2)| < K{\yi - 2/21 + \ti - t2|^/^). 

Theorem 10 For the stochastic differential equation (|lip . assume that (A4) is satisfied, andY^ of 
()12p is a sequence of random processes satisfying the local consistency condition. Then, the family 
of random variables {Y^{t) : t G [0, 1], /i > 0} are uniformly integrable, and Y^ =^ Y as /i — )■ 0. 

Proof: We divide the proof into several steps. 

1. First, we perform the following useful estimates. Set E^[-] := E[-|l^ = V^t!^ = t]^ ^^'^ 
z\t) := max{n >0:tt<t}, AM„^ = AY^^ - Ki^Yj^l 

Then 

E[supo<„<,.(t)_i I Er=o = E[supo<„<,.(i)_i I EtoKi^y!"] + AM/^I'] 



< KE 



sup I j;E,'^[AyA] 1 +Ke[ sup Ij^AA-//^ 



^0<n<zh{t)-l ' 



2i ~ ~ 
= /l + /2 
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We estimate the above two terms in the last hne separately. We have 

n „ 



sup ^Ef[AK,"] 



1 1 = E 

<e[ sup (Y,^^[^t^]{\h\{Y,\t^:) + 0{h))X by(LC2) 

<e[( ^\\^t\m{jt .t\) ^ o(K)) 

i=0 

<e[( j2 ^tH\bm\t'i)+o{h)) 



by Tower conditional expectation 



i=0 



< KE 



+ Kt^h^, 



i=0 
,z^(t)-l 



i=0 



i=0 



and 



n 

I2 =e[ sup I^AM^ 



i=0 
|2i 



-0<n<z''(t)-l 

< eH ^ Ai\/H j , by Doob's Maximal inequality 



i=0 



e[ J2 I^^/'I']' siiice E[(AM/^)'(AMj')] = 0,Vi/j 



i=0 



<i^E[ 5^ At1tv{aa'){Y^,t1) + {At'lf\b\'{Y!^,t1)]+Kht, by (LC2) (LC3), 



j=0 



where tr(A) denotes the trace of A. Therefore, we have 

E 



rh\2 



sup I V Ay/ 

-0<n<zhit)-l 

<i^(t + l)E[ A^/(|6|2 + tr(aa'))(y,^^,'^) 

i=0 



+ Kht. 



2. Applying inequality ([13]) . we have 



i=0 

<K\y\'' + KE[ Y ^tH\b\' + tT{aa')){Y:^,t1 



1=0 

t 



<K\y\^ + KE[ J {\b\^ + tT{aa'))(Y''{s),s)di 

<K + KE[f{Y^{s)fds\ 
Jo 



+ Kht 
+ K 
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Gronwall's inequality then yields 

E[|y^'(t)p] < K for all t G [0,T]. 
Thus we conclude {Y^{t) : < t < T, /i > 0} is uniformly integrable by Proposition [9l 



(14) 



3. We consider the tightness of {P/^ = P(y''') ^ : /i > 0}. It is enough to verify conditions 
imposed on Theorem 13.3.2 and the subsequent corollary given in [3]. 

(a) By Chebyshev's inequality and (fT^ 

lim limsupP{|y''(t)| > a} < lim lim sup ^E[|y''(t)|2] = 0. 



(b) For the purpose of characterization of tightness for discontinuous functions, we need to 
introduce some notions of modulus of continuity uj{Y^, 6) and ui'{Y^, 5) as the following. 
First, define uj{Y^,I) = sup^^^^^Gl \Y'^{ri) — Y'^{r2)\ for any subset I C [0, 1]. We also 
use TiS) to denote the collection of all (5-sparse partitions of [0, 1]. Then, we can define, 
for any 5 > 



and 



uj{y'',5)= sup uj{Y'',[t,t + 5]), 

0<t<l-5 



uj'{Y'',6)= inf max uj{Y'',[ti^i,ti)). 

{tjGr(5) i 



For the purpose of tightness of discontinuous functions, we also need to introduce mod- 
ified version of modulus of continuity. For an arbitrary e > 0, repeating the arguments 
of (fT3|) . we have 



lim lim sup F{uj'{Y'', 6) > e} 

< limlimsupP{a;(y'',25) > e} 

< lim lim sup ijElu'^CY'', 26)] 



< lim lim sup -^E 

< lim lim sup -^E 

< lim lim sup -^rE 



< lim lim sup K6 ( + h 

S-^O h-^Q 

As a result, {P,, = P(y'')~i : /i > 0} is tight. 



sup \J2^Y/'\ 

z'^(t)-l<n<z'^{t+25)-l 

n 

sup 5]A^,'^(|6|2 + tr(acT'))(y,^^, 

z'^{t)-l<n<z'^{t+25)-l 
t+2S 

Ki\Y''{s)\'^ + l)ds\ +Kh5 
0. 



+ Kh5 



1 
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4. Since is tight, for an arbitrary infinite sequence, there exists a subsequence that has a 
weak Hmit. For notational convenience, we denote this subsequence again by {Y^}, and its 
hmit by Y. Due to uniqueness of weak solution, it is enough to show that Y is the weak 
solution of (jlip . Since Y^ =^ Y, the uniform integrabihty and (A4) lead to 

hmE[Y''{t)]=E[Y{t)], and 
lim E[b{Y^{t),t)] =K[b{Y{t),t)]. 



Therefore, if we set 



we have 



M{t) :=Y{t)-Y{0)- [ b{Y{s),s)ds 

Jo 



E[M{t)] = limh-,oE[Y^{t) - Y^{0) - 



= lim e[ ^((yA - - / b{Yti,t)dt)\ - hm e[ / b{Y^\^^,s)d. 
Note that regarding the last term in the above, by linear growth of b 



(15) 



limE 



h 

zh(t) 



b{Y^\,^,s)ds\ 



< lim E 



'-Jt 



h 



K + K\Y^^^^^\ds < limE[At^,(,)(l + |y''(t)|)]. 



Since the Skorohod representation theorem (see [31 Page 70] ) allows us to work with a common 
probability space {Qi ,J^i,Fi) on which (At^^ , Y'^) converges to (0, Y) almost surely. Together 
with uniform integrabihty, it leads to 

Ihn E[At^,(^)(l + \Y\t)\)] = E[lim At^,(^)(l + |y''(t)|)] = 0. 



Therefore, we conclude that the second term on (jlSp satisfies 

ft 



limE 



'- Jt 



h 

h(t) 



b{YX,.,s)ds 



0. 



As for the first term of (|15p , using the Holder continuity in t and the tower property on local 
consistency 



iimE[^((yA-yAi)- r b{Y!L,,t)dt) 

1=1 -^^i-l 



_z^{t) 



hm e[ 5] {{Yt - Yt,) - 6(yA i,t)Att, + 0(|At, 



h |3/2 
i— 1 1 



)) 



i=l 



hm e[ {Atl,0(h^ + 0(|Atti|3/2)) 
ft— 5>0 L ^ — ^ 

«=1 



0. 
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Therefore, E[M(t)] = 0. In fact, one uses exactly the same procedure to show E[M(t)|J^5] = 
M(s) for any < s < t, and it concludes M{t) is a martingale process. 

Next, we use Yi and Yj^ to denote the ^th component of the vector process Y and Y^, and use 
{Yi,Ym){t) to denote the cross- variation of two real processes Yi and Y^ up to time t. Then, 

E\{Yi,Y^){t) - [\aiam){Y{s),s)ds\ 
Jo ^ 

= lim E|(y,\y^)(t)- / iaiam){Y\s),s)ds\ 

= \\Tamj2^Y|:,-Yl),_,){YX^-YX^-l)- I iaia^){Y\s),s)ds\=0. 

The last equality can be obtained similar to the above owing to E[M(f)] = 0, the local con- 
sistency, and regularity assumption (A4) on a. Therefore, the cross- variation {Yi,Ym){t) = 
/q (cr;crm)(?(s), s)ds in and in probability. This again implies that M{t) is a d-dimensional 
martingale process with its quadratic variation (M)(t) = f^{aa'){Y{s), s)ds in probability. 
Applying Levy's martingale characterization on the time-changed Brownian motion, there ex- 
ists a d-dimensional Brownian motion B{t) such that M{t) = fQa{Y{s),s)dB{s). Therefore, 
Y is the weak solution of (jlip . 

Summarizing all the above, it leads to the conclusion. □ 

3.2 Examples of Weak Convergence of Markov Chain Approximation 

The above construction of Markov chain approximation is based on the local consistency, which 
generalizes the local consistency in |14t Theorem 10.4.1] in the following senses. First, a and b may 
be unbounded but have linear growth. Therefore, the geometric Brownian motion is covered by 
weak convergence result of Theorem 1101 as an important application. Secondly, from the definition 
of local consistency in (9.4.2) of [13], Atjj has to be J-"^ measurable. However, the local consistency 
(LC1)-(LC3) allows us to have At^ to be only -F^+i measurable. 

In fact, locally consistent MC approximation is flexible for its various choices. For illustrations, we 
give several simple MC approximations for one dimensional process different from [T3], which can 
be easily verified to be locally consistent. 

Example 3 Euler approximation can be considered as a special case of MC approximations. Let 
{y^} be a Markov chain generated by 

1. y^f = y,t(j = 0. 
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2. Let the transition probability be Atjj; = h 

Fn'Vi = + 6(y„^ nh)h + a{Yj:,nh)VhNn, (16) 

where {A^n} are i.i.d. standard normal random variables. One can easily verify all local 
consistency, and Theorem 1101 implies . assuming (A4), the piecewise constant interpolation Y'^ 
of the form (jl2p is convergent to Y of (jlip . 

Example 4 The following MC approximation can be considered as an extension of binomial ap- 
proximation of Brownian motion, since it is adapted by approximation of W {t + h^) — W (t) ~ 
Bernoulli{\/2). Let d = di = 1. Let {Y^} be a Markov chain generated by 

1. Y^=y,t^^=Q. 

2. Let the transition probability is, with At^ = h? 

V^'Vi = Y:: + b{nh\ Yj:)h' ± a{nh\ = 1/2 (17) 

The above Markov chain is locally consistent since direct computation leads to 

1. E[sup„ At^] = h?, 

2. E[Ay^'|y„^ = y^tl = t]= n^t\yn = = • Ky.t), 

3. cov(Ay„^|y„'^ = = t)= mt'^,\Y^ = yA = t]- a\y,t). 

Therefore, assuming (A4), the piecewise constant interpolation Y'^ of the form (jl2p is convergent 
to y of ([II]) by Theorem Uni 

Example 5 This is a binomial tree type approximation of the diffusion term a{Y{t),t)dW{t) in 
(|lip . Let d = di = 1. Let {YJ^} be a Markov chain generated by 

1. Yq^ = x,4 = 0. 

2. Let AtJ^(y^,tfj) = ^, A /i^, and the transition probability is 

;y„\i = YJ: + 6(y„^^:J)A^:J(y„^^:J) ± (y„^^:J),a2(y„^^:J) > /i) = 1/2 

(18) 



y„V = >;" + K>;^^:^)A^::(y„^^:^) (y„^^::),a2(y„^^;i) < /ij = 1. 

Note that, this MC is also locally consistent since 
1. E[sup,At^] </i2, 
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2. n^Y,^\yn =y,t^n = t] = mti\Y:: = y,ti = t\- b{y,t), 

3. cov(Ay„'^|y„'^ = = t) = E[A4\YJ: = y,t';, = t]. {a\y,t) + 0{h)). 

Therefore, assuming (A4), the piecewise constant interpolation of the form ()12p is convergent 
to y of (dH) by Theorem Uni 



3.3 Can We Expect the Strong Convergence? 

In general, the answer is NO. To illustrate this point, we consider a special case of Euler approxima- 
tion of Example O Let {Q,J-,¥) be a probability space, on which J^t is filtration satisfying general 
conditions, and Wt of (|lip is a standard 1-d Brownian motion. We construct strong approximation 
of Euler-Maruyama's method by taking on iV„ of (jl6p by 

Wn+l - Wn 



Vh 



Under assumption (A4), the SDE (jlip has unique strong solution. Suppose Y'^ is a continuous 
interpolation of Euler approximation of Y'^ given by 

Y^{t) = Y^^ + b{nh, Y^^^){t - nh) + a{nh, - W{nh)), for t G [nh, nh + h). 

Then a classical result shows that see for example, |15^ Theorem 2.7.3], 



E 



sup \Y{t)-Y\t)\ 

0<t<T 



However, the above inequality fails for the piecewise constant interpolation of EM approximation 
{y'^}. Otherwise, we have following simple counter example. Consider EM approximation of Wt 
on [0, 1] by equal step size h = \/N . Then, we have 

,n — \ , 



E 



sup \W(t) - W{[Nt\/N)\ 

LO<t<l 



E 



sup 



sup 



\W{t)-W{- 



N 



■ l<n<N {n-l)/N<t<n/N 

Note that, W{t) = ^/NW{t/N) is a standard BM on time-scaled filtration. So one can reduce the 
above equality as 



E sup \W{t)-W{[Nt\/N)\ 
'-o<t<i 

where {A„} are i.i.d. random variables defined by 



1 



N 



E 



sup A„ 

l<n.<Af 



Kn = sup \W{t) - W{n - 1) 

n— l<t<n 



Since, A„'s are unbounded iid random variables, E 
shows that 

E 



supi<„<jv 



A^ 



goes to infinity as — t- oo. This 



sup \W{t) - W{[Nt]/N)\ > 0{N~^/^). 

0<i<l J 

In conclusion, one cannot expect more than weak convergence merely under local consistency. 
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4 Ramification and Further Remarks 



4.1 Application to Discretely Monitoring Barrier Option Underlying Stochastic 
Volatility Model 

We begin this section with an apphcation of Theorem [6] to the following stochastic volatility model; 
see [1]. Let W and B are two standard Brownian motions with correlation p in the filtered proba- 
bility space (fi, -F, P, . Suppose the stock price follows 

dX{t) = X{t)a{Y{t))dW{t) (19) 

with initial X(0) = x > 0, and volatility factor follows 

dY{t) = Y{t){fi{t)dt + b{t)dB{t)) (20) 

with initial 1^(0) = y > 0. We consider a discretely monitoring barrier option price given by ([3]). 

One can check that X and Y are unique nonnegative strong solution of SDEs ([19]) and (pOj) . if a 
satisfies polynomial growth, and /i and b are Holder- 1/2 continuous. In addition, we assume 

non-degeneracy of X, i.e., a{y) > for all y > 0- 

Note that it is possible to have cr(0) = under the above assumption. We examine the convergence 
lim/i = V when X^' is constructed by the Euler approximation. 

Note that the regularity condition (A2) is satisfied by [161 Proposition A.l] due to a{-) > 0. On 
the other hand, the payoff function 

g{xi, X2, • • • , Xfji, Xm+l, t) = {Xfn — I{ma.x{xi,...,Xm}>i} 

is only discontinuous at the points in the set {maxi X{i/m) = 1}. Thanks to the fact 

P{maxX(i/m) = 1} = 0, 

i 

it implies that g is almost surely continuous with respect to P. Another thing yet to be verified is 
the uniform integrability of {X^{t) : h,t}, since g is linear growth. Thanks to Theorem 1101 together 
with (A4), we have X^ X and the desired uniform integrability holds. Therefore, we reach the 
affirmative answer lim,, V'' = V. 

4.2 Further Remarks 

This work has been devoted to analyzing approximation to path-dependent functionals. We have 
obtained weak convergence of the approximation schemes, which provide a unified approach for 
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proving the convergence of numerical approximation of path-dependent functionals for a wide range 
of apphcations. 

This paper focused on diffusion models. For future work, it is worthwhile to examine systems 
driven by pure jump processes, jump diffusions, and systems with an additional factor process 
such as nowadays popular regime-switching processes. Systems with memory (time delays) form 
another class of important problems. Much work can also be devoted to numerical solutions of 
various stochastic differential equations, coordination of multi-agent systems, and many Monte 
Carlo optimization problems in which one needs to treat path-dependent functionals. 
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